Few-photon optical diode 



o : 

(N . 

3 : 

1—5 

m 



o . 
"—I . 

Oh- 

q ; 

O 

• I— I ■ 
Oh 



> 
oo 

rn 
O 

o 



X 



Dibyendu Roy* 

Department of Physics, University of California-San Diego, La Jolla, California 92093-0319 

We propose a novel scheme of realizing an optical diode at the few-photon level. The system 
consists of a one-dimensional waveguide coupled asymmetrically to a two-level system. The two or 
multi-photon transport in this system is strongly correlated. We derive exactly the single and two- 
photon current and show that the two-photon current is asymmetric for the asymmetric coupling. 
Thus the system serves as an optical diode which allows transmission of photons in one direction 
much more efficiently than the opposite. 

PACS numbers: : 03.65.Nk, 42.50.-p, 32.80.-t, 73.40.Ei 



A diode is an essential circuit element that allows uni- 
directional propagation of signal. Faraday rotator based 
on magneto-optic effect serves as an optical diode or op- 
tical isolator in most modern optical labs. There are 
also several proposed or demonstrated optical diodes re- 
lying on cascaded nonlinearities [1]. In this paper we 
propose a novel scheme to realize an optical diode at the 
few-photon level by making photons strongly interacting. 
Different schemes [2-6] have been explored to produce 
strong nonlinearities at the few-photon level. Recently 
a new approach has been suggested by Shen and Fan 
(SF) [7] to create photon-photon nonlinear interactions 
by coupling a bare two-level system (TLS) with a one- 
dimensional continuum for photons, such as a photonic 
crystal waveguide. Local interactions at the TLS induce 
strong correlations between photons. Based on a similar 
mechanism of realizing nonlinear interactions, an all op- 
tical single-photon transistor has been proposed by using 
surface plasmon modes of a conducting nano-wire [8] . 





SF consider a system of a one-dimensional waveg- 
uide with left and right-moving photons side-coupled to 
a TLS. The two-photon transport has been studied by 
them using a generalized Bethe-ansatz technique which 
properly takes into account the open boundary condi- 
tions of the problem [7] . It has been shown that the exact 
two-photon scattering eigenstates of this system include a 
two-photon bound state that passes through the TLS as a 
composite single particle. In this paper we study a similar 
electro-optical system that consists of a TLS directly cou- 
pled to a ID photonic waveguide (see top of Fig. 1). Also 
the coupling of the TLS with the waveguide is asymmet- 
ric. In a line-defect photonic waveguide the asymmetry 
in coupling may be realized by modulating the periodic- 
ity of photonic crystal at any side of the TLS locally. It is 
also possible to fabricate asymmetric coupling by work- 
ing with a system of superconducting transmission line 
resonators and a dc-SQUID based charge qubit [9, 10], 
or a system consisting of a bare two-level emitter inter- 
acting with the surface-plasmon modes of a conducting 
nanowire [8, 11, 12]. We construct the exact single and 
two-photon scattering states for our model by directly 
solving the single and two-particle Schrodinger equations. 



FIG. 1: (color online) (a) A schematic of the two- level system 
being directly coupled to photons in a one-dimensional wave- 
guide, (b) Unfolded model for the photons in the left (1) and 
right (2) side of the two-level system. 



Then we evaluate the single and two-particle current in 
the system. We find an interesting asymmetry in the 
two-photon current for asymmetric coupling between the 
TLS and the waveguide. This occurs due to a redistribu- 
tion of energy and momentum of photons after scattering 
from the TLS in the absence of spatial symmetry. Thus, 
our model acts as an optical diode which permits pho- 
tons to transmit in a direction more easily compared to 
the opposite direction. It is easier to achieve asymmetry 
in coupling for the direct-coupled TLS-waveguide system 
compared to the side-coupled TLS-waveguide. 

We consider a simple spin-boson type Hamiltonian 
to describe the direct-coupled TLS-waveguide system. 
The Hamiltonian, % = -I- J2i,k ^ikalk'^ik + 

Si.fe ^jfe(°^lfc + o,ik){.f^+ + o"-), where the first term rep- 
resents the free TLS with transition energy i7 (setting 
h—1); the second term corresponds to the free photons 
with energy Wj^, in the left (i = 1) and the right {i = 2) 
side of the TLS, and the last term describes the interac- 
tion between the TLS and the photons. Here aik (a|j.) 
is the annihilation (creation) operator for photons in the 
ith side of the TLS and Vik (V2fc) is the coupling be- 
tween the photon in the left (right) side of the TLS and 
the TLS. cr_ = ajoe (u+ = ajag) is the lowering (raising) 
ladder operator where aj(a|) is the creation operator of 
the ground (excited) state of the TLS. Then we unfold 
the model by performing standard manipulations for im- 
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purity models with fl being away from the cut-off of the 
linear dispersion (see bottom of Fig.l). Thus obtained 
the real-space Hamiltonian of the system as chiral ID 
field theories is given by 

H = —i ^ Vg / dx bl,{x)dxba{x) + ftaltte 

a=l,2 '' 

+ {Vib\{0)a_ + V2bl{Q)a_ + H.c.) , (1) 

where Vg [Vg] is the group velocity of the photons and 
b\{x) [b^ix)] is a photon creation operator at x in the 
left (right) side of the TLS. Vik {V2k) is replaced by a 
frequency independent Vi {V2). Here (T_|-|n, — ) = jn, +) 
and cr+|n, +) = 0, where \n, ±) represents the full system 
with n photons and the TLS in the excited (+) or ground 
(— ) state. Next we redefine the field operators in the two 
sides of the TLS, Ca{x) = ba{x). We assume here 

that Vg is real and a;— independent. Thus we find 

H = -i ^ / dx cl{x)dxCa{x) + ^lalue 

+ {Vic\{0)a- + V24(0)<T_ + H.c.) , (2) 

with Va = Val^Jv^- It shows that the renormalized 
coupling can be asymmetric for different group veloc- 
ities of photons at the two sides of the TLS even for 
Vi =¥2. We define the current operator in the system as 
/ = -ilH, Ni-N2]/2 where A''i {N2) is the total number 
of photons in the left (right) side of the TLS. Then we 
find / = i{Vic\{0)<7- - V2cl{0)a- - H.c.)/2. We want 
to find the expectation value of the current operator in 
the steady state. Thus we need to find the exact scat- 
tering eigenstates of the H. We here observe that under 
the transformation, ci{x) = {ViCe{x) + V2Co{x)) /V and 
C2{x) = {V2Ceix) - ViCoix))/V with V = sjvftv^ , the 
Hamiltonian in Eq.(2) breaks into two decoupled parts 
of even and odd field operators; i.e., "H = He + V-o 
where He = —i J dx cl{x)dxCe{x) + i}ala,, + V(cl{0)a- + 
<T+Ce(0)) and Ho = —i J dx cl{x)dxCo{x). In the trans- 
formed basis, the He contains the interaction term be- 
tween the even modes with the TLS while the Ho is just 
the kinetic energy of the noninteracting odd modes. The 
current operator in the even-odd basis is given by 



V 



^ct(0)a_ 



—^cl{0)a. 



H.c. 



(3) 



Single-photon dynamics: First we construct the exact 
single-photon scattering states of the full system and de- 
rive the corresponding single-photon current for a photon 
incoming from the left or the right side of the TLS. The 
single-photon eigenstates of the full system are 



/ 



photon is coming from which side of the TLS. is the 
excitation amplitude of the TLS. We find different coef- 
ficients in the eigenstates by solving the following equa- 
tions of motion obtained from the stationary single pho- 
ton Schrodinger equation. 



-idx9k{x) - kgk{x) + Vek6{x) = 
{n- k)ek + Vgkix)d{x) = 
-idxhk{x) - khk{x) = 



(5) 



where ^^(0) = [^^(0-) + gk{Q+)\/2 and gk{x) is a sin- 
gle photon plane wave for a; < 0. We find gk{x) = 



-x) + rkO{x) 



iV^/2), and 



hk{x) = 6'*=^, efc = V/{k -n + 

('k/c-k^ where 9{x) is the step func- 



tion. We define the single-photon scattering states as 
|1; fc) (|2; fc)) for those with an incoming photon of mo- 
mentum k from the left (right) side of the TLS. By choos- 
ing Ai = Vi/V, Bi = V2/V (Ai = V2/V, Bi = -Vi/V) 
in Eq.(4) we find the state \l;k) (|2; fc)). We then evalu- 
ate the steady state single photon current in |1; fc). 



J(l;fc) = (l;fc|/|l;fc) 



1 T/2T/2 

-Im CfeJ 



1/3 



TT (fc - fi)2 + (Li + La)' 



(6) 



with Fa = V^/2. As expected in the absence of any 
background phase shifts (which may arise if the TLS is 

embedded in a cavity) the single-photon current shows a 
Breit-Wigner-like (i.e., Lorentzian) line shape around the 
resonance k = Cl (see top two plots in Fig.2). For sym- 
metric coupling, i.e., Fi = r2 = F, a single photon at res- 
onance frequency is completely transmitted from one side 
of the TLS to the other side. Note here that the trans- 
mission and reflection coefficient for the direct-coupled 
model are similar to the reflection and transmission coef- 
ficient of the side-coupled model [7], a fact which also oc- 
curs for waveguide-resonator systems [13]. We also eval- 
uate (2;fc|/|2;fc) and we find |(l;fc|/|l;fc)| = |(2;fc|/|2;fc)| 
for arbitrary Vi,V2, i.e, the current is the same in mag- 
nitude for an incident photon coming from either side 
of the TLS. This even can be understood from the cur- 
rent expression in Eq.(6) which remains the same under 
exchange of Vi and V2 ■ 

Two-photon dynamics and optical diode: Next we 
construct the two-photon scattering eigenstates for this 
model. This is a nontrivial task. The two-photon incom- 
ing state for both the photons being incident from the 
left side of the TLS is given by 



dx 



{Ai[gkix)cl{x) + d{x)eka+] + 5i/ife(x)4(x)}|0, -),(4) 



dxidx2-^--^4>k{xi,X2)-^cl{xi)c{{x2)\0,-) , (7) 



with Ai and Bi being arbitrary constants to satisfy the 
nonequilibrium boundary conditions, i.e., the incident 



where ^k(a;i,a;2) = (e^'=i^i+^'==^= + eifeix2-Fife2xi) ^j^^ 
k = (fci,fc2)- We decompose the two-photon incoming 
state into ee, 00, and eo subspaces and determine the 
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full scattering eigenstates in the different subspaces sepa- 
rately. Note that to calculate the expectation value of the 
current operator in Eq.(3), we need to express the two- 
photon scattering eigenstate in the space of free photons 
as well as the TLS [14]. This is similar to the single pho- 
ton scattering states in (4) . Thus we differ here from the 
approach of SF [7] where the outgoing scattering states 
are expressed by a complete orthonormal basis spanning 
the free two-photon Hilbert space. 

The two or multi-photon transport in one-dimensional 
waveguides for the direct-coupled TLS model is strongly 
correlated. The local interaction at the TLS induces 
strong correlations between photons by preventing multi- 
ple occupancies of photons at the TLS. Thus the Hamil- 
tonian in Eq.l is equivalent to a bosonic version of the 
popular single-impurity Anderson model with an infinite 
repulsive interaction between photons at the TLS. When 
photons are incident upon the TLS, one of the photon 
gets absorbed and re-emitted by the TLS and then that 
interferes with the other photons in the waveguide. By 
the process of absorption and spontaneous emission the 
TLS mediates strong interactions between photons. The 
general two-photon scattering states are of the form: 

A2{g{xi, X2)^cl{xi)cl{x2) + e{xi)S{x2)cl{xi)(T+} 
+B2{j{xi:x2)cl{xi)cl{x2) + f{xi)6{x2)cl{xi)(T+} 

+C2 h{xi,X2)-^cl{xi)cl{x2) dxidx2\0,-) (8) 
v2 J 

where g{xi,X2) = g{x2-, Xi) and h{xi,X2) = h{x2,xi) due 
to Bose statistics of photons. e{x) (f{x)) is the probabil- 
ity amplitude of one photon in the e (o)-subspace while 
the TLS in the excited state. Here A2, B2 and C2 identify 
the boundary conditions for the incoming photons. Like 
the single photon case we here determine all the func- 
tions in Eq.(8) using the two-photon Schrodinger equa- 
tion with the total energy of the two photons E = ki + k2. 
Thus, we obtain the following linear equations. 



id. 



X2 - E]g{xi,X2) 
V 



e{xi)5{x2) + 5{xi)e{x2)\ 



V 



^ id^ - E + VL)e{x) + —[g{Q,x) 



id. 



X2 - E^i{xi\X2) + V6{xi)f{x2) 
idx - E + n^f{x) + Vj{{)]x) 



idx 



idx2 



E^h{xi,X2) = 



We solve the above equations with the initial condi- 
tions that the functions g{xi,X2), j{xi; X2), and h{xi,X2) 
in the region xi,X2 < satisfy the form of the cor- 
responding plane wave-functions given in Eq.(7). We 




FIG. 2: (color online) Plot of the single photon cur- 
rent [J(l;fci), 1(1; fc2)], the two-photon current change 
[5I{1; ki, fe), (5/(2; fci, fe)] with the transition energy Jl. Here 
fci = 0.3, k2 = 0.7, Vi = 0.45, and V2 = 0.15. 



also use, g{0,x) = 5(2;, 0) = [g(0+, a;) + .g(0— , a;)]/2, 
j(0,x) = [j{0+,x) + j{0—,x)]/2. In our solution we ex- 
plicitly separate out the two-photon bound state which 
arises in the ee-subspace. We find 



gixi,X2) 
e{x) 

j{xi;x2) 
/(^) 

h{xi,X2) 



^i(E-n+tV^/2)x,Q^^^ - X2)9{X2) + (1^2) 

^(5fci(a;)efe2 +5fc2(a;)efei) 

TT 

^ (.9fei {xi)hk2 (xa) + gk2 {xi)hk^ (xz)) 
■^(^ek^hk^ix) + ek^hk^x)^ 



(/)k(xi,X2) 



(9) 



The second terms in g{xi,X2) and e{x) are contribu- 
tions from the two-photon bound state. Both of these 
terms decay to zero with |a;i ~ X2\ or |a;| tending to infin- 
ity; this is an essential characteristic of the bound state. 
We emphasize that due to interactions at a localized re- 
gion two photons can exchange energy and momentum 
between themselves with the constraint of a fixed total 
energy and momentum. This redistribution in energy 
creates a strong interaction between photons. 

Now we calculate the expectation value of the current 
operator in the two-photon scattering state |1; /ci, k2) for 
both the incoming photons from the left of the TLS, i.e.. 
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FIG. 3: (color online) Plot of A/(fci,fc2) vs asymmetry in 
coupling X- Inset shows rectification TZ vs x- Here k\ = 0.3, 
fca = 0.7, Va = 0.15, and ^ = 0.5. 



A2 = V^/V^,B2 = ViV2/V\ and C2 = Vi /V^ 



(l;fci,fc2|/|l;A:i,fc2) 



27r2 1/3 



Im[efej + efej 



Im 



'fe2 fcl 



(10) 



Here C is the total length of the full system. The 
first term in Eq.(lO) is the contribution from two non- 
interacting photons and is similar to Eq.(6). The last 
part in Eq.(lO) is the change in two-photon current 
(call 5/(1; fcl, fc2)) due to photon-photon interactions and 
(5/(1; fcl, fc2) comes from the two-photon bound state part 
of the full wave-function. Also note that 5I{1] fci, fc2) is a 
factor C smaller than the noninteracting current; this sig- 
nifies that the probability of finding two photons near the 
TLS at any time is order of l/£ for a ID system of size C 
with two photons. C will disappear in the multi-photon- 
current in the system [15]. We can have W(l; fci, fc2) be- 
ing similar in magnitude of noninteracting current for an 
optimal two-photon interaction. This occurs when two 
photons are incident upon the TLS within a time de- 
lay of absorption and spontaneous emission by the TLS. 
We plot 5/(1; fcl, fc2) vs Vt for a fixed energy of incident 
photons in Fig. 2. We find from Fig. 2 that 5/(1; fci, fc2) 
changes sign near the single-photon resonance energy, 
i.e., = fcl or fc2. A positive 5/(1; fci, fc2) can be inter- 
preted as bunching of two photons after scattering where 
as a negative 5/(1; fci, fc2) is antibunching. 

We similarly derive the current for two photons be- 
ing incident from the right side of the TLS, i.e., A2 — 
Vi/V^,B2 = -ViV2/V^, and C2 = V^/V^. 



(2;fci,fc2|/|2;fci,fc2) = ^^^lu^lek. 



27r2 

V 

-Im 



2 * 



2 * 



(11) 



We see that the contribution from two noninteracting 



photons (the first part in Eq.(ll)) is similar in mag- 
nitude to that in Eq.(lO). But the change in two- 
photon current 5/(2; fci, fc2) (last part in Eq.(ll)) is dif- 
ferent in magnitude from 5/(1; fci, fc2) for Vi ^ V2. We 
plot —5/(2; fcl, fc2) versus Q. in Fig. 2. Hence we find 
that |(l;fci,fc2|/|l;fci,fc2)| ^ |(2; fci, fc2|/|2; fci, fc2)| for 
Vi 7^ V2^ i-e., the directly coupled ID waveguide-TLS 
system acts as an optical diode for two-photon trans- 
port. The mechanism behind diode like behavior in 
the present system is the redistribution of energy and 
momentum of photons after passing through the TLS. 
The asymmetry in the coupling generates asymmetry 
in the distribution of energy and momentum of pho- 
tons coming from the left or right of the TLS. In fact, 
|5/(l;fci,fc2)| : |5/(2;fci,fc2)| = /V^ is exactly same to 
the ratio of the bound state contributions for both the 
incoming photons from the left or right of the TLS. We 
remind that the redistribution of energy and momentum 
of scattered photons occurs via the two-photon bound 
states. A similar phenomenon for electrons has been 
studied recently [16]. Note that the present system acts 
as a diode in the fully quantum regime for a minimal 
two photons. It also differs from the classical spin-boson 
thermal rectifier [17] operating at a finite temperature 
bias. We here define A/(fci, fc2) — |(1; fci, fc2|/|l; fci, fc2)|~ 
1(2; fcl, fc2|/|2; fcl, fc2)|, and the dimensionless rectifica- 
tion, 7^ = 2A/(fci,fc2)/(|5/(l;fci,fc"2)| -f [5/(2; fcl, fc2)|). 
We also quantify the asymmetry in the coupling by 
X = 2{Vi - V2)l{Vi + V2). We plot A/(fci,fc2) and Tl 
with X for 8- fixed incident energy of two photons and a 
fixed Q. in Fig. 3. We use a constant V2 and change Vi in 
Fig. 3. We find that A/(fci, fc2) first increases with x ^^^d 
then falls after a maximum value. On the contrary, TZ 
increases monotonically with x ^-nd finally saturates at 
its maximum value 2. The multi-photon transport is also 
correlated in this system and it acts as an optical diode 
for many photons too. The technique employed here can 
be used to find current with three or more incoming pho- 
tons. 

The proposed optical diode may have potential appli- 
cations to build quantum circuits for optical quantum 
information processing and quantum computation [18]. 
However, the practical realization of the optical diode de- 
pends on achieving strong and asymmetric coupling be- 
tween the TLS and the photons in the waveguide. Strong 
light-matter interaction is the main focus of much re- 
search in quantum optics for quite sometime and has been 
demonstrated successfully in the recent past [9, 11, 19]. 
The asymmetry in coupling has not been much discussed 
before. The optical diode can be fabricated either with a 
TLS in a single waveguide or with a TLS being coupled 
to two different waveguides. It is technically challenging 
to fix a TLS at the right location in a single waveguide 
[9, 19, 20]. Again there may be high insertion loss for 
a TLS being coupled to two different waveguides. Re- 
cently an ultrahigh-Q coupled resonator optical waveg- 
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uides with low losses has been realized [21]. A strong 
direct coupling between an emitter and surface plasmons 
in a nanowire has also been achieved [11]. Hence it seems 
to be possible to reduce the insertion loss. In fact, the 
emitter-surface plasmons system is a promising candi- 
date for realizing asymmetry in coupling as the coupling 
between the emitter and surface plasmons is completely 
geometrical and depends on the diameter of the conduct- 
ing nanowire [12]. 

The discussions with S. Mookherjea, D. I. Schuster and 
A. Wallraff on the experimental realization of asymmet- 
ric coupling are gratefully acknowledged. The work has 
been funded by the DOE under Grant No. DE-FG02- 
05ER46204. 
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